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Abstract 

By considering the master equation of asymmetric exclusion process on a 
one- dimensional lattice, we obtain the most general boundary condition 
of the multi-species exclusion processes in which the number of particles 
is constant in time. This boundary condition introduces the various in- 
teractions to the particles, including ones which have been studied yet 
and the new ones. In these new models, the particles have simultaneously 
diffusion, the two-particle interactions AcAfj A^As, and the n-particle 
extended drop-push interaction. The constraints on reaction rates are ob- 
tained and in two-species case, they are solved to obtain a solvable model. 
The conditional probabilities of this model are calculated. 



1 Introduction 

One-dimensional asymmetric simple exclusion processes (ASEP) have been shown 
to be of physical interest in various problems in recent years. These problems 
are, for example: the kinetics of biopolymerization traffic models [3], 

polymers in random media, dynamical models of interface growth 0105 noisy 
Burgers equation [^j, and study of the shocks [Z1|H1- There are many review 
articles in these fields, see for example [Hj-fTK]. 

The totally ASEP has been solved exactly in ref.^HI- In this simple model, 
each lattice site is occupied by at most one particle and particles hop with 
rate one to their right-neighboring sites if they are not occupied. The model is 
completely specified by a master equation and a boundary condition, imposed 
on probabilities appear in the master equation. The coordinate Bethe ansatz 
has been used to show the factorization of the A^-particle scattering matrix to 
the two-particle matrices. 

By choosing other suitable boundary conditions, without changing the mas- 
ter equation, one may study the more complicated reaction-diffusion processes, 
even with long-range interaction. In ref. jl7|. the so-called drop-push model has 
been studied by this method. In this model the particle hops to the next right 
site even it is occupied. The particle hops to this site by pushing all the neigh- 
boring particles to their next right sites, with a rate depending on the number 

* alimohmd@ut . ac. ir 



1 



of right neighboring particles. The generaUzation of this model, by considering 
both the right and left hopping, has been done in ref-QH). Further generaliza- 
tion of the ASEP boundary condition have been proposed in As a 
comprehensive review on the ASEP, see |22| . 

All of the above studies are about the single-species systems. If one considers 
the model with more-than-one species, the situation becomes more complicated. 
The source of complexity is the abovementioned factorization of A^-particle scat- 
tering matrix, which in these cases restricts the two-particle S'-matrices satisfies 
in some kind of spectral Yang-Baxter (SYB) equation. In ref.jSHli all the solv- 
able two-species reaction-diffusion models, in which the number of particles is 
constant in time and the reaction rates are all equal, have been obtained. It 
was shown that there are 28 independent interactions (among 4096 possible 
types) which are solvable. One of these 28 models, in which the particles have 
exchange-interaction: A + B B + A, has been generalized to p species in 
ref.[2i]. In this model, besides diffusion to the right (with equal rates), particles 
interact through Aj + Ai ^ Ai + Aj with rate r.y . The spectrum of , to en- 
sure the solvability, has been obtained. Another multi-species reaction-diffusion 
model, which somehow relates to one considered in ,24,, has been discussed in 
P5| . In this model the particles have exchange-interaction with rates which are 
determined by the differences of their diffusion rates. 

The multi-species generalization of the models considered in |2Hj . has been 
considered in [201 ■ The processes are 



Aq0 -> $Aa with rate 1, 

• 



AaAp —> AjAs with rate . (1) 



Again the solvability restricts the reaction rates c"g. Some general remarks on 
the solution of SYB equation, and some special solutions of it have been given 
in ISSj. 

Recently, the generalization of the drop-push model to multi-species case has 
been considered in ref.|27]. The reactions are 



Aa9 %Aa with rate 1, 



AaAp^ -> ^A^As with rate . (2) 



The same drop-push reactions, but with n adjacent particles and with rates that 
are specified by a specific combinations of 's, also exist. Some comments have 
been given for solutions of SYB equation in The reactions (0) are called 
the extended drop-push reactions. 

In this paper we are going to study the most general multi-species model, i.e. 
the most general boundary condition, which all the previous mentioned models 
are the special cases of it. In its general form, the reactions are 



AJi 0A„ with rate 1, 

'^aAp — > AjAs with rate c"^, 

1)A^As with rate b"^ , 



(3) 



where the dots indicates the other drop-push reactions with n adjacent particles, 
in which in the meantime the types of the particles can also be changed. We 
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show that the reaction rates must satisfy some specific constraints, in order that 
we have a set of consistent evolution equations. The two-particle S'-matrices of 
this general multi-species model must also satisfy SYB equation, which in this 
case becomes more complicated. In the two-species case, we study a special 
class of reactions (jSJ in detail and show that in this class, there is only one 
solvable model, i.e. the solution of SYB equation is unique in this case. We 
study some physical properties of this unique model. 

The scheme of the paper is the following. In section 2, by using the law of 
conservation of probabilities, we obtain the most general boundary condition 
for a p-species reaction-diffusion model, in terms of two x matrices b and 
c. We obtain a constraint on the sum of the elements of each column of matrix 
b + c. We see that in special cases, this model is same as those have been 
studied previously. In section 3, we consider the matrix c as a diagonal matrix. 
It is shown that there are two cases. In the first case, the matrix b must be also 
diagonal, which the model reduces to the ordinary, i.e. single-species, drop-push 
model with variable rate. In the second case, c must be zero matrix, which the 
model becomes the one considered in [27| , that is the extended drop-push model. 
In section 4, we take c a non-diagonal matrix and obtain the necessary and 
sufficient conditions (constraints) for the matrices b and c, to have a consistent 
evolution equations. We see that the reactions are those indicated in eq.ljSJ. 
We find some classes of solutions of this set of constraints. In section 5, we 
investigate the coordinate Bethe ansatz solution for our problem and obtain the 
SYB equation, to ensure the solvability of the model. For p = 2, we show that for 
a specific class of parameters, which corresponds to a special exchange-reaction 
together with a special extended drop-push model, the SYB equation has a 
unique solution. Finally in section 6, the conditional two-particle probabilities 
are calculated for this unique model, and their large-time behaviour are studied. 



2 The general boundary condition 

Consider a p-species system with particles Ai, A2, ■ ■ ■ , Ap. The basic quantities 
that we are interested in are the probabilities Pai - aN i^i^ ' ' ' iXN',t) for finding 
at time t the particle of type ai at site xi, particle of type 02 at site X2, 
etc.. We take these functions to define probabilities only in the physical region 
xi < X2 < ... < xn- The most general master equation for an asymmetric 
exclusion process is 

d 

TTT^CKl ■ ■ - Qiv ('^1 : ' ' ' ; '^-/V 7 ^) -^Qi - ■■ajv('^l -^7 ' ' ' 7 "^N : ^" ' ' ' ^" -^ai ■■■ajv (-^l : ' ' ' )*^7V I?^) 

Ot 

-NPa-,...a„{xi,- ■ ■ ,XN;t). (4) 

This equation describes a collection of N particles drifting to the right with 
unit rate. It can be shown that if one decides to solve the master equation by 
coordinate Bethe ansatz method, one may only choose all the diffusion rates 
equal, which can be scaled to one [21]. So in eq.Q, we take all the diffusion 
rates equal to one. The master equation is only valid for 

Xi < Xj+i - 1, (5) 

since for Xi — x^+i — 1, there will be terms with Xi — Xi+i in the right-hand 
side of eq.Q), which is out of the physical region. One can, however, assume 
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that is correct for all the physical region Xi < ai^+i, and impose certain 
boundary conditions for Xi — x^+i. These boundary conditions determine the 
nature of the interactions between particles. Now the question is that what are 
the possible boundary conditions? To see this, we follow the same argument 
which has been given for single-species model in fSBJ . If one considers the master 
equations Q for two-particles probabilities, finds 

d 

X2 Xi<X2 X2 Xi<X2 



X2 Xi<X2 X X2 Xi<X2 

X X2 Xi<X2 

a; a: 

Now let us exclude the creation and annihilation processes, in other words we 
consider the processes in which the number of particles is constant in time. 
Therefore if we sum cq.(|H)l over ai and a2, the left-hand side becomes zero, 
so the right-hand side must be also zero. The only possible choice that leads 
the right-hand side of eq.® to zero, and consistent with more-than-two particle 
analogue of eq.® (see ref.[2HI) taking Pq^qj (a;, x\ t) as a linear combination of 
Pl3iP2 {x,x + \;t) and P/3i/32 (x — 1, x; <)'s. Therefore the most general boundary 
condition is 

P„,„, (x, x)=Y, bi\i\Pp,p, ix~l,x)+Y, cl\i\PM2 (x, X + 1). (7) 

(3 stands for (/3i/32)- These b and c matrices introduce interactions to particles. 
Inserting eq. (UJ in the right-hand side of eq. JHJ and sum over ai and a2 , results 



-EE ^"1"^ {x,x + i)+Y,Y.[T.(^+ ppiP^ix, x + i) = o. (8) 

X a X p \ a / 

Therefore the sum over the elements of each column of the matrix b + c must 
be equal to one, 

T.(b + c&2^1. (9) 

a 

Note that for annihilation process, the condition Q need not be satisfied |28|. 

If the matrices b and c are diagonal, then the boundary condition (0 does 
not induce reactions in which the types of the particles are changed. So the 
left-hand side of © is zero, without summing over ai and a2, and therefore 

- J2 P»io.2 {X, X + l)+ J2{b + c)l\ZPc.,c.2 {X, X + 1)=0. (10) 

X X 

So in diagonal case, we have, instead of condition 

b=l-c. (11) 
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The simple exclusion processes of QHl is an example of this model with p = 1 
and 6 = 0, the drop-push model with equal rate is p — I and c = case, and with 
non-equal rate, is an example of one-species case of eqs. ^ and l(TT|l |17j . The 
two-species model of [23 , special p-species of and general p-species reactions 
of |2ni are examples of eqs. and © with 6 = 0. Finally the extended drop- 
push model of [23 is an example with c = 0. Now we are going to consider the 
situations which have not been studied yet, that is the multi-species boundary 
conditions in which both b and c matrices present. 

3 Diagonal c 

As the first step of our study, let us consider the matrix c to be diagonal. As we 
will show in the next section, the non-diagonal elements of c are reaction rates 
and so they are non-negative, but the diagonal elements are not reaction rates 
and are defined through eq.®, which for diagonal c results 

« = i-Eev (12) 

a 

Now consider Pa^aii^i^ + !)■ Using eqs.Q, ijT)) and ((T^ . one finds 

P 

p 

The above evolution equation determines the following reactions as the source 
and sink of this model 

Aa% %Aa with rate 1, 

Ap,Ap,% %A^,A^, with rate (14) 

To check the consistency of our description, we next consider the 3-particle 
probability Paia2a3 {x—\,x,x+\). In this case, the boundary term Paia^ai {x — 
1, a;, a;) appears in the right-hand side of master equation. To calculate this term, 
we must use the boundary condition ITJ in two steps. The result is 

Paia2a3 — 1, X, x) = ^ ^ ^0203 ^ 1^2 "^Tl 72/33 i^X — 2^ X — \, x) -\- C^^^^-^aiftft — ^,X, x)] 

P 7 

"''''0203-^010203 (2^ ~ Ij a^j 2; + 1). (15) 

It is seen that the boundary term Pqiq2Q3 (a; — 1, x, x) is found as a linear com- 
bination of other boundary terms, i.e. ^01/32/33(2 — l,2,a;)'s. To avoid this 
problem, which in general does not lead to a consistent description of the inter- 
actions, we have two choices. One choice is to restrict ourselves to the cases in 
which the multiplication factors of the boundary terms in the right-hand side 
of cg. HlSI) arc zero. Then 

= (1 - E ^r^^rMit = 0- (16) 
7 
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So for non-zero reaction rates, i.e. ^ q, we have c"]^^ = (i.e. ^"172^ = 

1). In this case the model reduced to one studied in p7j, that is the extended 
drop-push model. 

The second choice is to consider the situations in which the boundary term 
in the right-hand side of ea. (|15|l is the same as one appeared in the left-hand 
side. In other words, we choose 

^Qi/32 "203 ~ ' "ia2a3"Q2"a3 ■ \^ ' ) 

This means that the matrix b must be also diagonal. Therefore we have 

Pal3{x,x) = rapPaffix - 1,X) + {1 - ra0)Pal3{x , X + 1) . (18) 

Considering Papix, x 4- 1), eqs.Q and (|18() show that the reactions are 

Aa^ ^Aa with rate 1, 

AaAf3% 9AaAf3 with rate Tq/j- (19) 

But if one considers Pa/3'y{x,x -I- l,a; + 2), or other n-adjacent particles proba- 
bilities, finds that the master equation Q), boundary condition lO, with b and 
c given in (|18(l can consistently describe the drop-push reactions only if Va/B is 
independent of a and /3. In this way we arrive at the one-parameter drop-push 
model, previously studied in TT . 

So in brief, considering c as a diagonal matrix, does not lead to any new 
model. 



4 Non-diagonal c 

Like the previous section, to find the reactions which exist in this case, we must 
first consider Paia2{x, x + 1). Using eqs.Q, and 0, we arrive at 

Pc.,c.Ax,X+l) = Pa,,c.Ax-l,X+l)+J2&,Pl3il32ix-l,x) 

+ E <'ti\Pp.P2 {x, x + 1)- 2P„,„, (x, X + 1) 
P 

= Po,,^,{x -l,x + \) + Y. b^\tPp.P.{^ -l,x) + J2 ci\tPp,0A^, x + 1) 

-(1 + E + E ^ + 1)' (20) 



in which we use 



13 fJ^a 

The evolution equation l|2()(l describes the following two-particle reactions: 



Ao,0 0Aa with rate 1, 

AaAjs —^ AjAg with rate c! 

AaA^9 iiA-^As with rate (22) 



Q/3 

7(5 ' 
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To find the more-than-two particles reactions, we must consider Pa-^...a^{x, x + 
1, • • • , a; + n — 1). In n — 3, we iirst need to know Pa-^a2a3{x — l,x,x). Using 
0, we have 

Paia2a3 ~ Ij ^7 ^) = ^ ] ^a^Jc^s ^')'l72/33 {x — 2, X — 1, x) + C^J^^P-yi^2l33 ("^ ~ "^i •^)] 

+ E cf^ia ^"1/32/33 - 1, X, X + 1). (23) 



Again we consider two cases. In the first case, we take the matrices b and c such 
that 

E 7i72t/32/33 _ XII Xf 2X133 (OA) 

Then eq.||2Sl) resuhs 

i,/3i/32 7,7/33 „/32/33 

(25) 

where o? = (ai, q;2, ■ ■ • )• Using (|25|l in computing Ps{x — 1, x, a; + 1), one can 
easily see that the resulting evolution equation gives different rates for same 
reaction. For example it gives the rate c^fa2 fo'" process /3i/32a3 ala2az^ 
which is consistent with second reaction of eq.(j22), and rate c^fa3/{^ ~ ''"a) for 
reaction ai/32/33 — ^ a\OLiOLj,. The only way to obtain the consistent result is 
taking = 0. Therefore we restrict ourselves to the matrices h and c with 
property 

E«^^f4=0- (26) 
Note that this constraint is equivalent to 

(l(g)6)(c«) 1) = 0, (27) 

in which h and c arep^xp^ matrices satisfying and 1 denotes thepxp identity 
matrix. Assuming and using eqs.Q), |3) and, (|^ . P^{x,x + + 2) 
becomes 

P5(a^, a: + 1,2^ + 2) = Pa{x-\,x + \,x^2) + Y.hif^^Pfi,p^o.3{x-\,x,x^2) 

13 

+ E Cola2J'ftfea3 (X, X + 1, X + 2) + ^ P^-(X - 1, X, X + 1) 

+ E ^t^Po^M {x,x + l,x + 2)-(l + Y. ^TpI + E Cftfc 
+ E + E Cfefe a; + 1, X + 2), (28) 

in which we use ea. lPT|l for d^l^l and cj^^2|. In above equation, 6^ stands for 

bl = Y.^it^lt3- (29) 
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The source terms of eq.(j2Hl), besides the reactions H22|l . are the following 3- 
particle drop-push reaction 

Ap^Af3^Af3^% ^ ^Aa^Aa^Aa^ with rate b'i. (30) 
The sink terms are also consistent with reactions H22(l and . provided 

T.^l-T.^r.%:k=T.^%\t:- (31) 

This is the last constraints that must be satisfied by the elements of matrix b. 
It can be shown that the more-than-three adjacent particles probabilities are 
consistent with following reactions 



Aa9 0v4a with rate 1, 

AaAfj ~* AjAs with rate c"g, 
^ao---^a„0 U^„---A^^ with rate (6„-i^„ • • • 6o,i)?o---7: ' (32) 

if the constraints lO, (|26|l and, (|31|l are satisfied. In ea. (|32|l . we use the following 
definition 

bk,k+i = 1 «) • • • (g) 1 (g) g)l g) • • • (g) 1. (33) 

k,k+l 

The important point is that we need not any further constraint. For example 
in describing the source and sink terms of the evolution equation of Ps{x,x + 
1, X + 2, X + 3), we encounter the constraint 

Eb$iT'^iXk = Y.f^r,:- (34) 

Now as ( see eq.ljSU) 

constraint (|^ reduces to J2(f3-y)^'ply^^'p2fy = E/3^/3i^2' ^^^^^ nothing but 
ea. (|31|l . This completes our investigation and the final processes are those in- 
dicated in lI221l- Of course, the question of solvability is open yet. 

As the elements of matrix b and the non-diagonal elements of matrix c are 
reaction rates, they must be non-negative, and this point restricts the allowed 
solutions of eqs.jnj, l|26|l and H31(l . Let us focus on two class of solutions. To 
be specific, we consider p = 2 case, but the arguments can be easily applied 
to arbitrary p. In p = 2, we may denote the two-particle states {ai,a2) as 
following: 

|1 (1, 1) , |2 (1,2) , |3 >= (2, 1) , |4 (2, 2). (36) 
solution of class 1: If b^^J^^ in ea. (|26|l is independent of (32, that is 

1,103 ^ /,2/33 ^ . . . ^ hP03 (Oy) 

then ea. l|^ reduces to 

E'=o;; = o. (38) 
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For p = 2, eq.(|27|) gives b^^ = b,^ ( for /Sg = 1) and 6^^ = 5^4 ( j^j. ^ 3), and 
ea. (|38|l for ai = 1, 2 specifies the elements of c as following 



/ 


-r 1 


r 2 



























-r 3 




V 








r 3 

C4 





(39) 



It can be easily shown that if we take the sum of elements of each column of b 
equal to one, then the constraints and H^lfl are also satisfied. The resulting 
model is a ten-parameter reaction-diffusion model. 

solution of class 2: If we take, for each /3, one or both of c^^J^ and b'^^as 
equal to zero, then ea. H26|l is satisfied. In p = 2, we can take c^^2^ — and 
^o!2Q3 = 0, which means that taking zero the second and fourth rows of c and 
first and second columns of 6. Then solving (pTT|l . results b2+b^ = 63^ + ^4* = 1 
and — X^i^j^- Finally considering jnj and noting the non-negativity of 

reaction rates, gives 



b = 



( 








\ \ 




( l-ci 


r 2 


r 3 










"2 






























, c = 






-r 3 





V 





i-bi 


1-&2* ) 




^ 








1 



(40) 
with condition 

{bi,b2\c,^)<l. (41) 

In p = 2, one can show that these two classes are the complete set of solutions 
of the constraint equations. Now we must seek the Bethe ansatz solution for 
the allowed set of parameters. 



5 The Bethe ansatz solution 

To solve the master equation (0J with boundary condition 0, we consider the 
following ansatz 

-fai , ■ ■ ■ ,ajv (^5 ^) ^ ■ ■ ,QAf (-'^) 1 (42) 

with 

^'(x) = ^A^e^'^fP)-''. (43) 

a 

^' is a tensor of rank N with components ipai,--- .qn (x), and the summation runs 
over the elements of the permutation group of N objects (SHllSni- Inserting 142(1 
in eqs.Q) and results, respectively, 

N 

EN^Y.i^-e-^P"), (44) 
fc=i 

and 

■ ■ ,Xk ^ x,Xk+i = X,- ■ ■ ) ^bkM+i'^i--- ,Xk^x~l,Xk+i=x,---) 

+ Cfc,fc+i*(- • • ,Xk ^ x,Xk+i ^ x + 1, - ■ ■ )(45) 
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Ck,k+i defined like bk.k+i in ea. H33|) . Tlic coefficients Aq. can be determined by 
inserting ea. (|43|l in (|45|l . which gives 

[l-e-'-(P^)6fe.fe+i-e'^(f'=+i)cfe,fc+i]A, + [l-e-'"(f'=+i)5fe,fc+i-e*'^(P'=)c,.fc+i]A,,, =0. 

(46) 

(Tfe is an element of permutation group which only interchanges pk and Pk+i- 

Uk ■■ {pi, ■ ■ ■ ,Pk,Pk+i, ■■■ ,PAf) -> (pi, • • ■ ,Pk+i,Pk, ■ ■ ■ ,Pn)- (47) 
Using H46|l . Acrafc is obtained from Ag- as following 

A^aak = Sk^k+lio-{Pk),(j{Pk+l))A.a, (48) 

where 

Sk,k+i{zi,Z2) = 1 «)•••«) 1 (g) S{zi,Z2)<»l «) • • • (g) 1. (49) 

k,k+l 

S{zi, Z2) is the following p^ x matrix 

Z2) = -(1 - Z^^b - Z1C)-1(1 - Z^H - Z2C), (50) 

and Zk = e'P*" . 

Until now, all the b and c matrices are acceptable. But we must note that 
the generators of the permutation group satisfy akUk+iUk — Uk+i'^kO'k+i, so 
one also needs 

In terms of 5- matrices, ea. H51|l becomes 

512(2:2, Z'i)S2'i{zi,Z'i)Si2{zi, Z2) = 523(^1, Z2)Si2{zi, Z3)5'23(^2, ^3). (52) 

Writing S'-matrix as the product of the permutation matrix 11 and a R matrix: 

Sk,k+1 ='■ ^k,k+lRk,k+l, (53) 

ea. (j52|l is transformed to 

-^23(22, 23)^13(2:1, 23)^12(21, 22) = i?12(2l, 22)-Rl3(2l, 23)^^23(22, 23). (54) 

This is the spectral Yang-Baxter equation |31|-|33j. 

Now the matrices b and c must be such that if one computes the two-particle 
S'-matrix H50|l by them, it satisfies the spectral Yang-Baxter equation H54(l . or 
equivalently ea. (|52|l . These equations are complicated and determining the ele- 
ments of b and c, in such a way that eq. (|52|l is satisfied, is not easy. 

For 6 = 0, the S'-matrix 1)50(1 becomes a binomial of degree one with respect to 
22 and therefore ((52(1 becomes a quadratic expression with respect to 23. Using 
this, one can transform the spectral Yang-Baxter equation 1(52(1 to a non-spectral 
matrix equation for matrix c, which its study is much easier than eg. 1(5211 |26| . 
The same is true for c = 0, since in this case eg. ((50(1 is linear in z^^ = e~^P^ 
|27( . When b and c are both different from zero, this procedure does not work 
because of the presence of 22 and 2^^ (or equally 21 and z^^) in both terms of 
(HB. 

For matrices b and c which are given in eg. l(40() . if one computes the cor- 
responding 5-matrix and writes eg. ((52(1 as RHS - LHS = 0, then one has 64 
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equations that must be solved for five variables bj'^, 62^7 ^ ^\ ^-nd . The so- 
lutions must be momentum-independent and non-negative. By obtaining the set 
of solutions of this system of equations by standard mathematical softwares, it 
is seen that all of the solutions are momentum-dependent except one. This only 
acceptable solution is: h'^ = 1, and other four parameters are zero. Therefore 
the solvable two-species reaction-diffusion model is defined through the following 
h and c matrices: 

/ \ 

10 


\ 1 / 

These matrices introduce the following reactions, with A = A\ and i? = A2, 



A% - 


IDA 


B% - 


%B 


AB - 


BA 


BA% - 


%AB 


BB% - 


-> %BB 


BAA% - 


%AAB 


BAB% - 


%ABB 


BBAfd - 


-> 9BAB 


BBB9 - 


^ 9BBB 



: (56) 

The dots indicates the more-than-three particle drop-push reactions which are 
specified by eg s . 1)321) and 1)33)1 . and all reactions occur with rate one. 



^ 1 
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(55) 



6 Two-particle conditional probabilities for re- 
actions ( 1561) 

The conditional probability P(c5, x; t|/3, y; 0) is the probability of finding par- 
ticles ai,a2, - ■ ■ at time t at sites xi, X2, ■ ■ ■ ^ respectively, if at i = we have 
particles /?2, ■ • • at sites yi, 2/2, ■ • • i respectively. In two-particle sector, it is 



Pab 
Pba 
\ Pbb J 



x;i|^,y;0) = / e-^=*e-^P''*(xi,x2)dpidp2 



(57) 



where by eqs.lgSJ, (|1H1), ISOI) and, (|S3), '^{xi,X2) is: 



^{Xl,X2) 



( " \ 

b 

c 

V d J 



e*P-^ + 5i2(zi,Z2) 



/ a \ 
b 

c 

V d ) 



2p.X 



(58) 
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in which p = (pi,P2), P = (P2,Pi), and 
by the initial condition, zi = ^ 



Z2 



stands for Ao-=i and specified 

\d ) 

'P^, and Si2{zi,Z2) is: 



512(21,^2) 



/ 

l-zi 

oof 



V 



22 









^2(1-Zl) , 
Zl(z2-1) / 



(59) 



After some calculations, one finds ( see for example refs.|21] and for more 
details) : 

P{A,A,^;t\A,A,y;0) = F,{t) + F2{t), 

P{A,B,x;t\A,B,y;0) = Fi{t), 

P(B,A,x;t|AS,y;0) = F^it), 

PiB,A,x;t\B,A,y;0) = Fi(i), 

PiA,B,x;t\B,A,y;0) = F^it), 

PiB,B,^;t\B,B,y;0) = Fi (i) + W , 



(60) 



-2t 



(xi - yi)! (x2 - ^2)!' 



and all other probabilities are zero. Fi{t)^s are: 
F^it) ^ * 

i^4(t) - 
J^5(0 = 



'2t 



{xi-y2 + iy- {xi-y2) 
{x2 - yi - 1)! (a;i - y2)!' 

p2-l/l fXi-y2 + l 

{x2 - yi)! (xi - 2/2 + 1)!' 

f^x2-yi ^x2—yi — l 



E 

fe=0 



(2:2 -yi + k)l'' 



(a;2-yi)! (a;2- 2/1-1)! 



E 

k=l 



fXi-y^+k 



(61) 



Note that beginning with initial state {A, B), the system can go to state {B, A) 
by third reaction of (|56|l and again goes back to {A, B) by fourth reaction, etc.. 
As the rates of reactions are equal, we expect that at large time, the probabilities 
of finding {A, B) and {B, A) are equal, if they are summed over all accessible 
sites. In fact, if one calculates the difference of these two probabilities, finds: 

00 X2 — 1 

E E [P{A,B,^;t\A,B,y;0)^P{B,A,x;t\A,B,y;0)] 



X2=y2 xi=yi 



2 J2 Im{2t) + Io{2t) 



(62) 



in which In{x) is the n-th order Bessel function of the first kind. Now at x ^ cxd, 
we have 

In{x) ^ (63) 



'2lTX 



12 



so 

M 

lim D{t) = 0, (64) 

t^oo v47rt 

where M is the number of /„(2i)s in ea. (|62|l . The final result confirms our 
expectation. The same result can be proved when the initial state is (B, A). 
Acknowledgement: I would like to thank N. Ahmadi for useful discussions. 
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